Fluids in extreme confinement 
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For extremely confined fluids with two-dimensional density n in slit geometry of accessible width 
L, we prove that in the limit L — > the lateral and transversal degrees of freedom decouple, and 
the latter become ideal-gas-like. For small wall separation the transverse degrees of freedom can 
be integrated out and renormalize the interaction potential. We identify nL 2 as hidden smallness 
parameter of the confinement problem and evaluate the effective two-body potential analytically, 
which allows calculating the leading correction to the free energy exactly. Explicitly, we map a fluid 
of hard spheres in extreme confinement onto a 2d-fluid of disks with an effective hard-core diameter 
and a soft boundary layer. Two-dimensional phase transitions are robust and the transition point 
experiences a shift 0(nL 2 ). 



Confined fluids are intermediate between fluids in three 
and lower dimensions. The confinement strongly influ- 
ences their physical behavior [1] like structural [2] and 
dynamical [3] properties, and in particular the phase be- 
havior and phase transitions [4, 5] or the glass transi- 
tion [6]. Therefore, confined fluids have attracted a lot 
of attention during the last three decades. One of the 
widely discussed features is the influence of the restricted 
geometry on the critical behavior (see Ref. [4, 7-10] and 
references therein). For a colloid-polymer mixture in a 
slit geometry with walls separated by five colloid diame- 
ters the critical exponents are already very close to those 
of the 2d- liquid [7]. Reduction of the spatial dimension 
from three to two replaces in the Kosterlitz-Thouless- 
Halperin- Nelson- Young (KTHNY) theory [11] a first- 
order phase transition by a two-stage continuous transi- 
tion from a fluid to a hexatic phase and then to a 'solid' 
with long-range orientational, but algebraically decaying 
translational order [2, 12]. 

Most studies consider a slit geometry with two par- 
allel, hard plates separated by a distance L + a with a 
single component fluid of either hard spheres with di- 
ameter a or point particles (a = 0). We will consider 
the case of extremely confined fluids where only a mono- 
layer fits between the walls. Most investigations address- 
ing this regime used computer simulations (see Ref. [13— 
15] and references therein), density functional theory 
(DFT)(see Ref. [16-19]), virial expansion, free volume 
theory, effective diameter theory [14] , and integral equa- 
tions [20, 21]. For instance, the density profile at the 
center of a neutral hard-sphere fluid between two paral- 
lel neutral hard walls (HSHW) has been calculated ex- 
actly for Lja — > [19, 20]. How the 3d-functional for 
the excess free energy reduces to the corresponding 2d- 
functional was investigated for a HSHW for < L < a 
in Refs. [16, 17]. The approximate analytical approach 
and the computer simulations [13, 14] as well as exper- 
iments [22, 23] of a HSHW reveal that the transition 
of the 2d-system from the fluid to the 'solid' triangular 
phase and from the latter to a buckling phase, a stable 



bilayer phase first observed and theoretically explained 
in refs. [22, 24], exists up to Lja ~ 0.6. For Lja > 0.6 
new phases appear which do not have an analogue in the 
'2d- world'. That the hexatic phase (probably not found 
in Refs. [13, 14] because the system is too small) exists 
for i/o"Lj < 0.15 has been shown recently for a Lennard- 
Jones (LJ) liquid of particles with a 'diameter' olj [15]. 

Intuitively it is obvious that a fluid confined between 
two plates approaches a 2d-fluid as L — > 0. But one of the 
interesting questions is: How does such a fluid converge 
to a 2d-fluid if the effective distance L becomes smaller 
and smaller? Or vice versa: If a 2d-fluid, e.g. undergoes 
an equilibrium phase transition, what is the range of L 
such that the transversal degrees of freedom (d.o.f.) do 
not affect the properties of that transition? To study 
these questions is the main motivation of the present 
work. Surprisingly, analytically exact results can be de- 
rived, which is a rare situation for strongly interacting 
many-particle systems. We will show that for L — > the 
lateral d.o.f. decouple from the transversal ones, where 
the latter behave as an ideal gas in the external wall po- 
tential. Additionally we will calculate the exact leading 
order correction due to their coupling. This allows us to 
determine the leading L-dependence of thermodynamical 
quantities. 

We investigate a fluid of N identical particles with lat- 
eral and transversal d.o.f. Xj = (r i; z^), i = 1, 2, . . . , N. 
The fluid is confined between two plates at z = ±(L + 
<t)/2 parallel to the x-y plane, and the area of a plate is 
denoted by A. The particles mutually interact via a pair 
potential only 

= (i) 

i<j 

where we abbreviate Xij := \xi — Xj\. To illustrate 
our approach we consider a pure hard-core repulsion 
with core diameter a and ignore additional particle- 
wall interactions. A generalization to smooth pair 
and particle-wall interactions is straightforward. As 
usual the configurational partition function reads Z = 
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FIG. 1. Illustration of the cluster expansion: The hard core 
with diameter ctl is shown in gray and the dashed circle with 
diameter a marks the range above which the excluded vol- 
ume interactions vanish. 2d-fluid of hard disks with diameter 
ctl (a), a 2-cluster (b), two irreducible 3-clusters and an irre- 
ducible 4-cluster (c) immersed in the hard disk fluid. 



Tr (exp [~ f3V ({xi})]) with the configurational integral 
Tr (') = / [n,=id 3 xj (•)■ Then F ex = -k B T\n(Z/V N ) 
is the excess free energy [25] with respect to a three- 
dimensional ideal gas of accessible volume V = AL. 

Let us outline the strategy of our approach. First, it 
will be shown that the configurational part p({xt}) — 
Z~ x exp [— /3V({a!i})] of the canonical ensemble factor- 
izes for L — > into a transversal and lateral distribution. 
In a second step a cluster expansion with respect to a 2d 
reference fluid is developed, which allows us to eliminate 
the transversal d.o.f. and to obtain an effective potential 
V c s({fi}; L), that adds to the bare potential of the ref- 
erence fluid. Third, wc employ Vrfj({rj}; L) to calculate 
the leading correction to the free energy. 

The key observation is that by pure geometrical rea- 
sons the lateral coordinates of two neighboring spheres 
cannot come closer than a L = v 'a 2 — L 2 for small L 
as follows directly by Pythagoras' theorem. Hence the 
effective interaction after tracing out the transversal de- 
grees of freedom is concentrated on a thin shell of area 
7r(cr 2 — <j\) = ttL 2 in addition to the bare two-dimensional 
hard-core repulsion of diameter a L , to be denoted by 
VF({r,}; L) = J2i<j W(rij;L). To avoid cumbersome no- 
tation, we mostly suppress the explicit dependence on L 
in the following for the interaction potentials. Then the 
distribution function (as a measure) factorizes to lead- 



ing order p{{$i}) = p±({z i })p\\({f f i })[l + 0(L 2 )] into 
the distribution function of the lateral and transversal 
d.o.f.. Here the transversal distribution reduces to a 
one-dimensional ideal gas, p±({zi}) = 1/Z±, and triv- 
ial partition function Z± = L , whereas the lateral 
d.o.f. correspond to a two-dimensional hard-disk fluid 
P||({^}) = ex P( — PW({ri}))/Z\\. In particular, the free 
energy simplifies to F = F id + fL + 0{L 2 ) where F id 
is the three-dimension ideal gas contribution due to the 
kinetic energy, and fJL = — fc^T hx(Z\\ /A N ) is the excess 
free energy of a hard-disk fluid of diameter <j L . 

Next, we elaborate the leading correction to the fac- 
torized ensemble. Let us introduce the cluster functions 



= e(r£ + (Zj - z,) 2 - a 2 ) - 6(rf 



with Q(x) the Heaviside function, and := |r*, — fj\. 
Note, —1 < fij < and its support is restricted to 
°l < < a. Then the identity exp[— /3V({xi})] = 
exp[— /^^({rj})] ni<i(l + fij) ai l° ws us to perform a 
cluster expansion. Let us emphasize that the subse- 
quent procedure can be directly generalized to the case 
of additional smooth pair and particle-wall interactions 
or point particles by suitable choice of the reference 
potential ^({ri}) and cluster functions fij. For con- 
venience we abbreviate pairs i =/= j by a — (ij) and 
enumerate them. Then, we define the effective poten- 
tial exp(-/3V r eff ) = (l\ a {l + f a ))± by averaging over the 
transversal d.o.f. . With x = i^\ a {l + f a ))± — 1 and us- 
ing the series expansion ln(l + x) = J2^ =1 (—l) k+1 x k / k, 
one finds 



-f3V cS =J2 



■ E KiW»± - ifa)±(fp)±] 

a<f3 



(2) 



The first line contains precisely the terms 
(—l) k+1 J2 a [(fa)±} k /k and adds up to the exact 

(2) 

effective two-body potential V^ ff , the subsequent 
term contains the first contribution to the three-body 
interaction, 

N 

V cS = -k B Tj2 M 1 + (/a)x) + Yl ^(k-clusters). (3) 

a k=3 

Note that successive cluster contributions are additive 

(k) 

and involve fc-body interactions V^ s , see Fig. 1 for il- 
lustration of the various clusters. Keeping only the two- 
cluster term is equivalent to the approximation (J} Q (1 + 
fa))± ~ l\ a ( l + ifa}±)- For hard spheres {f a )± can be 
evaluated explicitly 



Vg>(r ir ,L) = -2k B Thx{l ^/{a 2 r 2 3 )/L 2 ), (4) 
for ctl < Tij < a and zero otherwise. The total pair po- 

(2,} 

tential Vtotal(»") = W(r) + V^ s '(r), represented in Fig. 2, 
smoothly interpolates between the hard-core repulsion 



3 



of disks with diameter <jl and the force-free region for 
r > a. The additional effective interaction diverges log- 
arithmically for r I Oh and approaches zero as a square 
root yj (a 2 — r 2 )/L 2 for r t o - . 

The effective potential can be used for the calcula- 
tion of ensemble averages. Consider an observable X = 
X({ri}) which depends only on the lateral coordinates. 
Then its configurational average yields 



PV total (r) 



3 PA F{Xe -ev e{f) 



(5) 




where AF := F - F id - Fl = -fc B Tln(exp(-/3Kff))|| 
is the shift of the free energy induced by the coupling 
between transversal and lateral d.o.f., and (-)|| indicates 
a configurational average with respect to pn. Hence the 
knowledge of the effective potential allows calculating all 
structural quantities of the lateral d.o.f. in the confined 
system. In particular, using the effective two-body po- 
tential yields averages that are correct including to order 
0(nL 2 ). 

The shift of the free energy AF can now be evaluated 
explicitly to leading order. Since the support of is 
concentrated to a tiny shell of width ~ L 2 , the corrections 
with respect to the reference system become small and 
induce a hierarchy of contributions of decreasing weight. 
Then the average (exp(— PV e f[))ii, c.f. Eq. (2), evaluates 
to a power series in L 2 . Abbreviating y = ((J| Q (1 + 
/ )) x )|| - 1, yields —(3AF = ^ =1 (-l) fe +VA- Since 
y = 0(L 2 ) , we find for the leading correction of the free 
energy 



/3AF = ]T((/ Q K)||+0(L 4 ). 



(6) 



The leading correction arises from the effective pair po- 

(2) 

tential . In the thermodynamic limit one infers 

AF/N = J g(r)[e-^S^ -l]d 2 r + 0(nL 2 ) 2 , 

(7) 

where g(r) is the radial pair distribution function of the 
hard-disk reference fluid. Since the effective potential 
acts only on the thin layer, g(r) can be replaced by its 
contact value g(<J^) and the remaining integral be per- 
formed with the result 



AF/N = —Trk B TnL 2 g{a A 



0(nL 



2\2 



(8) 



The preceeding relation uncovers the smallness param- 
eter nL 2 of the confinement problem, which is one of 
our principal results. Let us emphasize, that this co- 
incides only formally with a virial low-density expan- 
sion; the quality of our approximation does not arise due 
to a dilute system, rather by the strong confinement. 



FIG. 2. The total two-body potential as function of distance 
r. The gray region represents the excluded volume region 
corresponding to the closest lateral distance ox two spheres 
can assume in confinement. 



The contact value is connected to the excess surface 
tension S ox = — (dF ex / '&A)t,n of the two-dimensional 
reference system via the virial equation E cx /n/csT = 
■nna\g(a^) /2, similar to the three-dimensional case [25]. 

Let us discuss some consequences of these results. 
For instance, the force per area exerted on the plates 
p(T,L,n = N/A) = —A~ 1 (dF/dL)T,A,N follows to 



P 



nk B T 
L 



1 



G 



imL 2 g(a + ) + 0{nL 2 f 



(9) 



where the leading term is of purely entropic origin and 
arises from the ideal gas term in the transversal direction. 
The excess free energy f] x of the reference hard-disk sys- 
tem does not contribute and the coupling of the lateral 
d.o.f. to the transversal ones is evaluated to leading or- 
der. Similarly, the surface tension 



nk B T 



= 1- 



5nnL x 
12 



g(a+)+0(nL 2 ) 2 , (10) 



consists of the corresponding surface tension of the ref- 
erence hard-disk system and the corrections due to the 
coupling. Since the effective diameter of the reference 
system crj, depends also on the plate distance, the cor- 
rection due to the shift of the hard-disk surface tension 
is of the same order as the correction due to the coupling 
to the transversal d.o.f.. 

As an application of our findings, assume that the ref- 
erence 2d-fluid undergoes a phase transition at a two- 
dimensional packing fraction <p*\ As shown above, the 
leading correction to the bare potential of the 2d refer- 
ence fluid is of the order L 2 . Consequently, for small 
L we have ^ 2d) (i) = yf d) {L = 0)[1 + 0(L 2 )}. Note, 
that this remains true if besides the hard core inter- 
actions smooth pair and particle-wall interactions are 
added, as well as for point particles. The result for a 
HSHW [13, 14, 23] allows us to quantify this behav- 
ior. The 3d and 2d packing fractions are related by 
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FIG. 3. Phase diagram of a HSHW taken from Refs. [13, 14]. 
Symbols represent the MC data points, the dashed lines are 
guides for the eye and the thin solid lines represent L((p( 3d ')/a 
from Eq. (11). 1 A and b denote the triangular and buckling 
phase. Inset: Phase transition lines ip( 2d ^ as function of L/a. 

(p( 3d )(L) = (2/Z)^ 2d \L)/{l + L/a). Then the above 
discussion leads to 

V i 3d) (L) ~ 4 3d) (L = 0)/(l + L/a). (11) 

Therefore, to leading order the L-dependence of the 
phase transition line ipi 3 ^ (L) arises only from the trivial 
factor 1/(1 + L/a). The corrections are 0(L 2 ) and orig- 
inate again from the coupling of transversal to lateral 
d.o.f.. Figure 3 shows part of the phase diagram ob- 
tained in Ref. [13, 14] including our leading order result, 
Eq. (11). The figure demonstrates that the freezing and 
melting line between fluid and triangular phase are well 
described by our analytic prediction up to L/a < 0.3. 
The freezing phase boundary between triangular and 
buckling phase follows Eq. (11) even up to L/a < 0.5. 
Since the buckling phase develops a transversal struc- 
ture with increasing L/a, our cluster expansion cannot 
be applied, because macroscopic clusters will be involved. 
The range of L/a for which the 2d-behavior dominates 
becomes even more evident using the 2d packing fraction 
(c.f. the inset of Fig. 3). This inset clearly demonstrates 
the shallow rise of the curvature of the transition line 
with increasing L/a and their vanishing slope at L = 0, 
as predicted by our analytical result. Additionally, the 
MCT data hint that the slope of the melting line of the 
b-phase is nonzero at L = 0. 

To summarize and conclude, we have proven that the 
lateral and transversal d.o.f. of an extremely confined 
fluid in slit geometry decouple if the effective width L be- 
comes much smaller than the average lateral particle dis- 
tance 1/y/n. Since the transversal d.o.f. approach ideal 
gas behavior for L — > 0, the nontrivial thermodynamic 
properties are completely determined by the correspond- 
ing 2d-fluid. The leading correction to the free energy 



due to the residual coupling has been calculated exactly, 
thereby identifying nL 2 as smallness parameter of the 
problem. The phase behavior in extremely small slits 
is close to the underlying 2d-fluid and we conclude that 
phase transitions are robust. Let us emphasize again, 
that our approach is valid for a densely packed strongly 
interacting system in strong contrast to an ordinary virial 
expansion of a dilute gas. Beyond thermodynamics all 
structural properties can be evaluated correctly in next- 
to-leading order in an effective two-dimensional ensemble 
where the two-body interaction consists of a hard-disk 
repulsion with reduced diameter and a thin smooth re- 
pulsive layer. Due to this mapping a simulation of the 
hard-sphere fluid in extreme confinement could be re- 
placed by a simulation of an effective 2d-fluid. Similar 
conclusions apply for the construction of functionals in 
density functional theory (see Ref. [26]). 

Phase transition lines are analytic in the vicinity of 
L = 0. Within the radius of convergence of our cluster 
expansion no phase transition can occur, however it is 
unclear if this radius in general signals a morphological 
transition. In particular, the analytic property implies 
the existence of a critical width L c , in case that the 2d- 
fluid does not exhibit a single first order transition as 
found for a fluid of hard disks [27, 28] (see also the dis- 
cussion in Ref. [29]). The two transition lines emerging 
from both 2d transition points have to join into a first 
order line at L c , independent on whether the KTHNY 
scenario or that of Ref. [30] holds. Whether L c is finite, 
as found for the 3-state Potts model [31] which has a con- 
tinuous transition in d = 2 and a first order transition in 
d = 3 is not yet obvious. Recent experiments for colloidal 
films may support a finite critical thickness L c [32]. 

Let us conclude by stating that it is straightforward 
to include a smooth wall potential and a smooth part on 
top of the hard-core repulsion. Then the reference en- 
semble still factorizes, the transversal partition function 
being a product of single-particle contributions. Corre- 
spondingly the force on the plates acquires an additional 
contribution from the wall potential. Our framework can 
be applied also to point particles with smooth interac- 
tion potentials and suitable wall potentials. Then the 
effective potential is proportional to the mean-square dis- 
placement ((Az) 2 )j_ of the transversal d.o.f. which is of 
the order L 2 and the free energy shift is again of the or- 
der nL 2 . In particular, from our analysis one can design 
specific particle-wall interactions minimizing ((Az) 2 )i to 
stabilize a two-dimensional phase behavior, e.g. the hex- 
atic phase investigated in Ref. [15] for attractive walls. 

The applicability of our approach and its results are 
universal for any confined c£-dimensional fluid where the 
effective width L of one of the spatial extensions con- 
verges to zero. The corresponding d-dependent small- 
ness parameter is given by the dimensionless quantity 
nL d_1 , where n is now the number density of the d— 1- 
dimensional fluid. 
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